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Within a theoretical model based on the Boltzmann equa- 
tion, we analyze in detail the structure of the unusual peak 
recently observed in the relaxation time in Cu. In partic- 
ular, we discuss the role of Auger electrons in the electron 
dynamics and its dependence on the d-hole lifetime, the op- 
tical transition matrix elements and the laser pulse duration. 
We find that the Auger contribution to the distribution is 
very sensitive to both the d-hole lifetime and the laser 
pulse duration ti and can be expressed as a monotonic func- 
tion of ti/th- We have found that for a given Th, the Auger 
contribution is significantly smaller for a short pulse dura- 
tion than for a longer one. We show that the relaxation time 
at the peak depends linearly on the d-hole lifetime, but in- 
terestingly not on the amount of Auger electrons generated. 
We provide a simple expression for the relaxation time of ex- 
cited electrons which shows that its shape can be understood 
by a phase space argument and its amplitude is governed by 
the d-hole lifetime. We also find that the height of the peak 
depends on both the ratio of the optical transition matrix ele- 



ments R = \Md^s 



and the laser pulse duration. 



Assuming a reasonable value for the ratio, namely R = 2, 
and a d-hole lifetime of th = 35 fs, we obtain for the calcu- 
lated height of the peak Ar t h = 14 fs, in fair agreement with 
At eip « 17 fs measured for polycrystalline Cu. 



I. INTRODUCTION 

Two-photon photoemission (2PPE) has been applied 
intensively to study the electron-electron interaction in 
alkali, noble and transition metals-oH The experimen- 
tal results are. often compared with the Fermi-liquid 
theory (FLT)j3 but such a comparison is only valid if 
the influence of secondary-electron generation and trans- 
port of excited carriers is negligible. Indeed, a series 
of 2PPE experiments involving photoexcitation of elec- 
trons from the d band in Cu have shown an urwsual non- 
monotonic behavior of the relaxation time.LTc2l The re- 
laxation time shows a broad peak at energies just above 
the d-band threshold which shifts linearly with photon 
energy. Both the non-monotonic behavior and the de- 
pendence on photon energy are incompatible with FLT 
description or with prQerit calculations based on Green's- 
function formalism.EJTj There is currently a controversy 
in the literature on the origin of the-ueak and especially 
on the role of secondary electrons JjHj One of the key 
questions concerns the amount of Auger electrons con- 
tributing to the hot-electron distribution in the region 
of the peak. On the basis of simple model calculations, 
Petck et al. have estimated that the Auger contribu- 
tion is less than 10%£3 while Knocsel, Hotzel, and Wolf 



could reproduce their correlation measurement by assum- 
ing that the hot electron population above the d-band 
threshold was entirely generated via the Auger process.0 
Note that the experiments were performed using different 
pulse durations (12 and 70 fs, resp.). 

In a previous work, we have developed a theory for 
the dynamics of excited electrons in metals, which in- 
cludes both the effects of secondary electrons and trans- 
port. We used a Boltzmann- type equation in the random- 
Is approximation for the calculation of the distribution 
function of excited electrons. A detailed description is 
given in Ref. The calculations have shown a peak in 
the relaxation time at the right energy reflecting the d- 
band threshold and a linear shift with photon frequency, 
in agreement with experiments. However, the height of 
the peak Ar t h ~ 10 fs calculated using Th =n 9 fs was 
smaller than the observed one Ar exp « 17 fs.tJ Further- 
more, a simplified procedure was used to extract the re- 
laxation time which also led to an overestimate of the 
peak height. Indeed for simplicity, the relaxation time 
was extracted by fitting the distribution of excited elec- 
trons, f(t), which was very sensitive to the influence of 
secondary electrons causing a deviation from a simple 
exponential decay. To make our results more reliable, in 
this work we use the same fitting procedure as the one 
used by experimentalists, i.e. we fit / 2PPE (Ai), given by 
the convolution of the probe laser profile with the distri- 
bution of excited electrons. 

In this paper, we will analyze in detail the structure 
and height of the peak and the role of secondary elec- 
trons, especially focussing on the Auger contribution and 
the d-hole lifetime. In ocder to shed light on the contribu- 
tion of Auger electrons ,EHl3 we have started by analyzing 
their weight in the total hot-electron distribution. For 
this purpose we have performed calculations for differ- 
ent hole lifetimes (t/j) and laser pulse durations (17) and 
found that their relative contribution can be expressed 
in a simple way as a product of a function of 17/77, and 
a phase space function which is energy-dependent. In 
addition, we show that the height of the peak in the re- 
laxation time depends mainly on the hole lifetime and 
surprisingly not on the amount of Auger electrons. It 
will be shown that the relaxation time for energies above 
the d-band threshold scales linearly with Th and its struc- 
ture is determined by the same phase space function as 
mentioned above. Finally, we study how the height of the 
peak depends on the optical transition matrix elements 
and on the laser pulse duration. 
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II. RESULTS AND DISCUSSION 

As said in the introduction, the composition of the 
hot-electron distribution in terms of different contribu- 
tions (f(E) = J2i fi(E), i = P,A, OS for primary, Auger, 
other secondary electrons) is an open question. We de- 
termine their respective weight and compare the distri- 
bution of hot electrons for different cases in Fig. [j]. We 
especially focus our attention on the variation as func- 
tion of the hole lifetime at the top of the d band, Th- 
The hole lifetime plays an important role, since it de- 
termines the generation rate of Auger electrons which 
may strongly influence the dynamics of excited electrons 
above the Fermi energy. Note that in our model, the 
hole lifetime is the inverse of the Auger scattering rate, 
thus we cannot vary independently the hole lifetime and 
the amount of Auger electrons generated. One should 
mention that experimentally, lower bounds of 24 fs and 
26 fs have been reported-for the hole lifetime at the top 
of the d band in CuOtJ Since there is up to now no 
clear consensus on the value of Th, in our calculation we 
will consider it as a free parameter. In absence of a d 
band, the electron-hole symmetry would hold, and our 
calculation gives then = 9 fs for the hole lifetime at 
E — Ep — —2.2 cV. Due to the localized nature of the 
d-wave functions, the d-hole lifetime is expected to be sig- 
nificantly larger than the s-hole lifetime, which will hence 
be considered as a reasonable lower bound. We start with 
the case of a laser pulse of photon energy hv = 3.3 eV and 
duration n = 70 fs (FWHM) as used in Ref. § Later, 
we will discuss the effect of changing the laser pulse du- 
ration. 

In Fig. |l|(a) we show results of the distribution of 
hot electrons for different cases: with and without sec- 
ondary electrons and for different hole lifetimes. We 
observe a peak in the primary electron distribution at 
E — Ep = 1.1 eV which is due to the sharp feature in the 
(i-band density of states; at this energy, primary electrons 
are dominating. To simplify the discussion, in Fig. |l|(b), 
we plot fi/ fp for i = A, OS. In the region of the peak 
in /, the primary electron contribution is dominating, 
which appears now as a dip in fi/ fp. Below the peak, 
there is a continous rise of both the Auger and other sec- 
ondary electron contribution. On the other hand, just 
above the peak, there is a sharp increase of the weight of 
the secondary electrons, largely dominated by the Auger 
contribution, with a maximum at E — Ep = 1.3 eV. We 
observe that a change of the hole lifetime leads to a signif- 
icant variation of the Auger contribution. For example, 
for t/j = 35 fs, /a//p = 2.2, whereas for shorter lifetime 
Th = 9 fs, = 3.6. In our model, a shorter hole life- 

time corresponds to a larger Auger scattering rate and 
therefore to a larger amount of Auger electrons. Note 
that /a//p has the same shape for any given Th, only 
the relative amplitude is changing. This will be analyzed 
further when discussing the effect of the laser pulse dura- 
tion. Furthermore, we find that for E — Ep > 2.5 eV, the 



contribution of secondary electrons is always negligible. 

In Fig. ^, we show the relaxation time as function of en- 
ergy for different cases (see Fig. |l|). The relaxation time 
is determined by a fit to the 2PPE signal as discussed in 
connection with Fig. ||. It will be interesting to see how 
the peak in the relaxation time depends on both f aI fp 
and on the hole lifetime. First, we observe that when 
considering only primary processes {t{E) = Tp(E)), the 
relaxation time is a monontonic function, as it can be 
expected on the basis of FLT. However, when including 
other secondary processes (but no Auger process), we im- 
mediately observe some significant changes for E — Ep < 
2.3 eV and some new structure at intermediate energy. 
When comparing with Fig. 0(b), the changes in t(E) due 
to the presence of other secondary electrons coincide with 
the variation of fos I fp discussed previously. Indeed, 
at E — Ep = 1.1 eV (dip in fos/ fp), T is almost un- 
changed. Then a small peak appears at E — Ep = 1.3 eV 
(region where the secondary contribution is the highest), 
and for E - E F > 2.5 eV, t(E) = t p (E) (negligible 
secondary electron contribution). This picture remains 
qualitatively unchanged when including the Auger pro- 
cess. 

Let us now comment on the quantitative aspect. First 
of all, a clear and well-defined peak is observable when 
Th gets sufficiently large (typically Th > 17 fs). It is 
also interesting to note that for small Th (large Auger 
scattering rate) , the effect is very weak, in spite of the fact 
that the Auger contribution to the electron distribution 
is large [see Fig. |j](b)]. This is a clear indication that 
the amplitude of the peak in r is mainly controlled by 
the hole lifetime. This is indeed illustrated in the inset 
of Fig. Ij, where r at E — Ep = 1.3 eV is plotted as 
function of the hole lifetime. We see clearly that t scales 
linearly with Th- 

It is important to mention that we have found that this 
linear scaling behavior is also valid for other energies. 
This suggests that t(E) could be expressed as a sum 
of two terms, tq(E) (including primary and other sec- 
ondary processes), and a term proportional to Th (Auger 
contribution), where the coefficient is a phase space fac- 
tor depending on the DOS only. To be more explicit, 
we suggest that the relaxation time can be written as 
t(E) — t (E) (x T h F(E), where F(E) is a phase space 
factor defined as the ratio between a) the available phase 
space for Auger scattering and b) the phase space for the 
optical excitation. 

To test the validity of this expression, in Fig. ||, we 
have plotted in the same figure [t(E) — To(E)]/Th as 
function of E — Ep and the calculated F(E). Clearly, 
for E — Ep > 1.3 eV, one can see that the values for 
[t(E) — To(E)]/Th lie almost on the same curve, which 
depends only on the energy. A deviation is observed for 
smaller energy, where Auger electrons are no longer the 
dominant contribution to the hot-electron distribution 
[see Fig. [|(b)]. However, the position of the dip is always 
found at 1.1 eV. Within some deviations at high and low 
energy, we observe a fair agreement between F{E) and 
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[t(E) — to(-E)]/t/j. Note that a peak appears in F(E) at 
E-E F » 1.7 eV, whereas in [r(.E) — to(.E)]/t/i, one finds 
only a shoulder at this position. This can be attributed 
to the fact that by definition F(E) contains the explicit 
details of the DOS, whereas the relaxation time is not so 
sensitive to the details of the DOS. 

To summarize this discussion, the important result is 
that the amplitude of the broad peak in t(E) is controlled 
by 77, and its shape by the phase-space factor F(E) (in 
random-k approximation, it depends only on the DOS). 
When comparing Fig. || and Fig. [l|, we find that Ja/Ip 
is also proportional to F(E), which will be commented 
later. 

As was shown in Fig. ||, the shape of the peak in the re- 
laxation time is determined by the function F(E) which 
contains the information on the optical transition matrix 
elements. Therefore, in Fig. ^ we analyze how sensitively 
the relaxation time t(E) depends on the optical transi- 
tion matrix elements. Till now, we have presented data 
for the case of equal matrix elements, Md-> S p — M sp ^ sp . 
However, on the basis of the different degree of localiza- 
tion of the d and sp-type wave functions, it is expected 
that Md^ S p should be larger than M sp ^ sp . For example, 
in the calculation of the dielectric functicp-,for Ag, a ratio 
\M d ^ p \ 2 /\M p ^ s \ 2 = 2.21 was estimated.^ 

In Fig. |[ we have plotted t(E) for different values 
of the ratio R = \Md^ sp \ 2 /\M sp ^ sp \ 2 . If we define the 
height of the peak as At = r(1.3 eV) — r(1.2 eV), we find 
that when changing R from 1 to 2, At varies significantly 
from 9 fs to 14 fs (increase by 60%). However for a fur- 
ther increase of R from 2 to 4, the change is much weaker 
(14 fs to 17 fs, increase by only 30%). Assuming that for 
Cu, R is of the same order of magnitude as for Ag, one 
can get an estimate of the hole lifetime required to get 
fair agreement with the observed height of the peak for 
polycrystalline Cu, AT exp sa 17 fs.B We have found that 
for R = 2, t% = 35 fs gives Ar t h = 14 fs. Note that our 
results obtained in the random-k approximation are most 
suitable for comparison with experiments performed on 
polycrystalline material. Measurements on single crys- 
tals for different orientations have provided values rang- 
ing from A:zu_ = 15 to 40 fs, which is the same order of 
magnitude. LrE3 It is interesting to note that 77, = 35 fs 
coincides with the value, used by Knoesel, Hotzel, and 
Wolf in their simulations 

Note that our conclusions are not restricted to the case 
of Cu only, but hold also in the case of other noble and 
transition metals. Indeed, let us discuss the case of Ag, 
Au and the 3d transition metals Fe, Co and Ni. In the. 
case of Ag, no peak is observed in the relaxation time,l!l 
due to the fact that the d band threshold is approxima- 
tively at 4 eV below the Fermi level and no d electrons are 
excited for the photon energies widely used in 2PPE ex- 
periments (hv < 4 eV). The case of Au is more intriguing, 
since the location of the <i-band threshold is very similar 
to the one in Cu. Indeed, a very weak structure in the 
relaxation time was also observed for Au.E2l The fact that 
the peak in Au is so weak can be attributed to the small 



degree of localization of the of the d hole. It is expected 
that a 5d hole in Au should be less localized and there- 
fore have a larger Auger scattering rate than a 3d hole in 
Cu. Thus one expects r^ u < r^ u and thus in accordance 
with our result that the peak is governed by t/,, a less 
pronounced structure in Au than in Cu. Another possi- 
ble explanation is the fact that the peak in the <i-band 
DOS in Au is much less pronounced than in Cu (again 
due to the less localized d electrons in Au), leading to a 
less sharp separation between excitations from the d and 
sp bands and therefore to a weaker peak. Let us now 
comment on the fact that in the ferromagnetic transition 
metals Fe, Co and Ni, no peak is observed. U tQ There are 
two main reasons for the absence of this feature, (i) A 
threshold for d — > sp transitions is not observable due to 
the band structure. The majority d band extends to en- 
ergies above Ep and a d — > sp threshold does not exist. 
The upper minority d band edge is located very closely 
to Ep and the threshold occurs only at the highest exci- 
tation energies where secondary (Auger) electron contri- 
butions are not important, (ii) The d holes have a large 
available phase space for decay within the d band which 
leads to a very short ci-hole lifetime. 

All the results presented up to now were for 77 = 70 fs. 
However, experiments were performed withJaser pulse 
durations ranging from 70 fs down to 12 fs .LJ E£j Hence, 
it is interesting to calculate how t{E) depends on 77. 
First, let us show that the relative weight of the Auger 
contribution to the hot electron distribution is indeed 
strongly sensitive to the laser pulse duration. In Fig. 
we have plotted Ja/Ip as a function of l/r/j for laser 
pulses varying from 12 to 120 fs duration. The first 
observation is that when reducing the laser pulse du- 
ration for a given r/j, /a/Zp is strongly reduced. For 
instance, for Th = 35 fs, Ja/Jp goes from 3 to 0.5 (6 
times smaller) when reducing the laser pulse from 120 
to 12 fs. This might already explain the origin of the 
controversy about the amount of the Auger electron con- 
tribution mentioned before. It is clear that in experi- 
ments with very short laser pulses (12 fs in Ref. |lj), the 
Auger electron contribution is much smaller than in ex- 
periments with longer laser pulses (70 fs in Ref. ||). At 
this level, one can conclude that the laser pulse dura- 
tion T] is a new relevant time scale. However, we have 
found that the only relevant parameter for the variation 
of /a//p is in fact 71/77,. In the inset of Fig. |5[ we have 
plotted /a//p as function of 71/77, and we see that all 
the data points lie to a very good approximation on the 
same curve, i.e. fA/fp(n,Th) oc 5(77/77^. In this figure, 
the data are plotted for E — Ep = 1.3 cV, however, it has 
to be stressed that this scaling is also found for higher 
energies. The ratio /a//p can be reduced to an expres- 
sion of the form f A /fp{E,T h T h ) oc g(n/Th)F(E), where 
F(E) is again the phase factor discussed in the previous 
section. 

Since the relaxation time is determined from 
7 2PPE (At) which contains explicitly the information 
about the laser pulse, let us now analyze how r is af- 
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fected by changing the pulse duration. In Fig. |[ we 
show the dependence of the relaxation time t(E) on the 
pulse duration for durations varying from 120 to 12 fs. 
We observe a significant effect in the region of the peak. 
Interestingly, the peak gets higher for the shorter pulses, 
although the amount of Auger electrons is much smaller, 
see Fig. |5[ Note that the dependence on the pulse dura- 
tion disappears for high energies, where the secondary- 
electron contribution is negligible. In order to under- 
stand the sensitivity to the pulse duration, in the inset of 
Fig. ||, we present the calculated / 2PPE and the fit by an 
exponential model function used to extract the relaxation 
time. We see that the agreement is very good for pulses 
of 120 and 70 fs duration, while there is a clear deviation 
from the single-exponential fit for shorter pulses of 40 and 
12 fs. Indeed, for the 12 fs pulse, there is a delayed rise 
in / 2PPE which is due to the delayed generation of Auger 
electrons. Such a delayed rise for verv_short laser pulses 
was recently observed experimentally.EEl As an additional 
remark, the observed deviations underline the difficulties 
in extracting a relaxation time in the limit of very short 
pulses and indicate that maybe an improved definition of 
the relaxation time would be appropriate in this case. 

III. CONCLUSION 

To conclude, this work presents a detailed analysis of 
the role of secondary electrons in 2PPE experiments on 
Cu. The d-hole lifetime plays a crucial role. Our con- 
clusions are more general and can also be applied to 
other noble and transition metals as discussed. We have 
found that in the case of Cu, the secondary-electron dis- 
tribution is dominated by the Auger contribution (for 
Th ~ 35 fs) for longer laser pulses (77 > 40 fs), but not for 
a very short laser pulse of r/ = 12 fs. For a given Th, the 
Auger contribution is much larger for a longer laser pulse 
duration than for a shorter one. The parameter which 
controls the variation of Ja/Jp is tj/t/j. Concerning the 
structure in the relaxation time, we have shown that the 
relaxation time at the peak depends linearly on Th, but 
surprisingly not on the amount of secondary electrons 
generated. The shape depends on a phase-space factor 
(in random-k approximation, on the DOS). We provided 
an expression for the relaxation time t(E) as a sum of 
two terms: the first for primary (and other secondary) 
electrons and the second for the Auger contribution. We 
have also found that the height of the peak depends sen- 
sitively on the optical transition matrix element ratio, 
R = \Md^sp\ 2 l\M sp -, sp \ 2 , and on the laser pulse dura- 
tion. For a value of R — 2 and Th — 35 fs, the calculated 
height of the peak is At^ = 14 fs in fair agreement with a 
measurement on polycrystalline Cu giving Ar exp w 17 fs. 
Note, preliminary results indicate that transport changes 
the magnitude of the relaxation time, but not the struc- 
ture and height of the peak. The influence of transport 
will be studied in detail in a separate publication] 2 .!! In 



view of the importance of the <i-hole lifetime, it would be 
highly desirable to perform further experiments as well 
as theoretical calculations on the d-hole lifetime in both 
Cu and Au. 
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FIG. 1. a) Distribution f(E) of hot electrons at t=0 for 
different cases (1: only primary, 2: primary and other sec- 
ondary, 3-5: primary, other secondary and Auger electrons 
for different values of the hole lifetime th)- b) Ratio of distri- 
bution of Auger (/a) and other secondary electrons {fos) to 
distribution of primary electrons (fp). The pulse duration is 
ti — 70 fs and the photon energy hv = 3.3 eV. 
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FIG. 2. Dependence of excited-electron relaxation time 
(determined by a fit to the 2PPE signal) on hole lifetime th 
as a function of energy. The inset shows the relaxation time 
at the peak position, E — Ef = 1.3 eV, as a function of Th- 
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FIG. 3. Ratio [t(E) — ro(E)]/r h discussed in the text, 
where to(E) includes only effects of primary and other sec- 
ondary electrons (curve 2 in Fig. ^). F(E) is the phase space 
factor discussed in the text. 



FIG. 5. Ratio of Auger to primary electron contribution, 
Ia/Jp at E — Ep = 1.3 eV, as a function of the inverse hole 
lifetime 1/t;, for different laser pulse durations r;. In the inset, 
we show /a//p as a function of ti/th- The line is a guide to 
the eye. 




FIG. 4. Dependence of the relaxation time on the ratio of 
optical transition matrix elements R = \Md->s P \ 2 /\M sp ^ sp \ 2 . 
The laser pulse duration is t; = 70 fs and the hole lifetime is 
n — 35 fs. 
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FIG. 6. Relaxation time for different laser pulse durations 
ti. The hole lifetime used is = 35 fs. In the inset, we show 
the calculated two-photon photoemission intensity J 2PPE a t 
E — Ef = 1.3 eV (solid line) together with the fit assuming ex- 
ponential decay from which the relaxation time is determined 
(dashed line). 
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